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Abstract 

We present a theory of the current-voltage characteristics of a magnetic domain wall between 
two highly spin-polarized materials, which takes into account the effect of the electrical bias on the 
spin-flip probability of an electron crossing the wall. We show that increasing the voltage reduces 
the spin-flip rate, and is therefore equivalent to reducing the width of the domain wall. As an 
application, we show that this effect widens the temperature window in which the operation of a 
unipolar spin diode is nearly ideal. 

PACS numbers: 



The discovery of the giant magnetoresistance effect |l| and the rapid growth in the number 
of its industrial applications have raised the hope that a similar breakthrough, perhaps of 
even broader consequence, may result from the combination of established semiconductor 
technologies with a precise control of the spin degree of freedom. As part of a growing effort 
in what has been called "semiconductor spintronics" [2|, |3( several new spin-based devices 
have been designed and discussed during the past few years: we mention, for example, the 

Datta-Das A \ spin transistors, the bipolar spin transistors of Zutic and Das Sarma |5| and 

□ 

Flatte et al. p and, lastly, the unipolar spin diode and transistor of Flatte and Vignale [7|, 



. All these devices, while still largely theoretical, are actively pursued in the lab, since 
they might eventually prove useful for computer operation such as nonvolatile memory and 
reprogrammable logic. 

At the heart of the above-mentioned devices is a magnetic junction (or magnetic domain 
wall), i.e., a region of inhomogenous magnetization connecting two regions of different homo- 
geneous magnetizations. In this paper we extend the conventional theory of spin transport 
across such a junction to include the effect of the electric field in the inhomogeneous region 
between two highly spin-polarized materials. Our work is motivated, in part, by recent 
insights on the role of electric field on the efficiency of spin injection across a magnetic 
interface l9( and, more specifically, by the recent discussion of the unipolar spin-diode in 
Refs. jj], |S|] . A simple model for this device is two ferromagnetic conducting slabs, denoted 
F\ and F2, with oppositely aligned magnetizations, connected by a domain wall of width d. 
The direction of the exchange field B(x) within the domain wall rotates linearly through an 
angle n in the z-x plane, i.e., 

B{x) = B [cos9(x)x + sin9(x)z], (1) 

where —n/2 < 9 < 7r/2 and < x < d. We distinguish between the component of the 
current due to "up-spin" electrons, J-f, and that due to "down-spin" electrons, J^, and 
accordingly define the charge current J q = Jj + J| and spin current J s = J-j- — Jj_, where 
"up" points in the positive x direction. If the domain wall is sufficiently sharp (i.e., more 
precisely, if d is much smaller than , h , where m* is the effective mass of the electrons 
and A is the magnitude of the exchange splitting between the up- and down-spin bands) 
then the spin of an electron crossing the junction is essentially conserved. Under these 
conditions a unipolar device (where the charge carriers on both sides have the same polarity) 



is analogous to a classical p-n diode, with up- and down-spins corresponding to electrons and 



holes, and t 
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ie oppositely aligned magnetic regions playing the role of the p-type and n-type 



materials [7(. A bipolar device (where the charge carriers on different sides have opposite 
polarity) can also be analyzed in this context under conditions of forward bias. A key 
assumption, particularly in the analysis of the unipolar spin diode, is that the applied bias 
voltage drops almost entirely across the junction, whose resistance is therefore supposed 
to be much higher than that of the rest of the structure. Indeed, recent experimental 
work has confirmed that highly resistive and well localized domain walls can be realized at 
nanoconstrictions in GaAs HU Coherent spin transport aeross hi ghly resistive vertical 

nnn 

tunnel junctions [12J, [13J, [14j may also be analyzed based on models such as we present here. 
An important deviation from ideality, namely the possible occurrence of spin-flip processes 
in the junction, was examined in detail in Ref |2J. Such spin-flip processes are responsible 
for the appearance of a lower critical temperature below which minority-spin injection is no 
longer operative and direct tunneling between the majority-spin bands perverts the operation 
of the diode. However, the analysis of Ref. [q] did not account for the electric field that is 
present in the domain wall region when an external bias is applied. From the high-resistivity 
assumption we know that this field is significant, and from the work of Yu and Flatte [9| we 
know that even a modest electric field, in a semiconductor, can have a large and favorable 
effect on the efficiency of minority-spin injection. These considerations motivate us to refine 
the analysis of |a] to include the effect of the electric field on the spin-flip rate. The outcome 
of the improved analysis is both interesting and reassuring: on one hand, it shows that the 
electric field greatly favors minority-spin injection, thus widening the temperature window 
in which the spin-diode exhibits an "ideal" behavior; on the other hand it confirms the 
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essential validity of the original treatment of Ref. 

We now review some of the essential aspects of the analysis from which the results above 
are obtained. In pursuing the natural analogy between p-n diodes and unipolar spin diodes, 
a number of assumptions are required, which closely correspond to those introduced by 
Shockley for an ideal diode [l5j |: (1) within the diode, the voltage drop occurs mainly across 
the domain wall junction, (2) the Boltzmann approximation for transport is applicable, (3) 
the minority carrier density is small compared to that of majority carriers, and (4) there 



is no "recombination current" in the domain wall. It has been argued in 



that (1) holds 



if the domain width is sufficiently small (in the sense specified above). Additionally (2) 



holds if the voltage is not excessively large, (3) if the spin splitting is large compared to the 
temperature, and (4) if the spin coherence time is much larger than the time required to 
traverse the domain wall. 

With these assumptions in mind, we can begin a reconstruction of the I—V characteristics 
by considering the action of a single electron incident on the domain wall. There are four 
possibilities [from the four possible combinations of reflection (r) or transmission (t), with 
spin flipped from its original orientation (sf) or not flipped (nf)], the probabilities of which 
will be denoted: r s f, r n f, t s f, t n f. Throughout our analysis, we will consider this set of 
coefficients to be the controlling quantities in the behavior of the spin diode, as they form 
the basis for all subsequent calculations. When a voltage V is applied to the diode, we 
can think of the regions -F\ and F 2 as two majority spin reservoirs of opposite alignment at 
quasi-chemical potentials /xi = and /12 = eV, respectively, which, it has been observed, 
are not appreciably altered by the presence of current. Then the majority- and minority- 
spin currents in these regions, due to electrons with energies in the range (E, E + dE), are 
described component-wise by (see ref. ja|): 

j n {E) = -(l-r nf (E))f n (E)+t sf (E)fo(E) 

ME) = r sf (E)f n (E)+t nf (E)fo(E) 

j^(E) = (1 - r nf (E))fo(E) - t sf (E)f n (E) 

j 2l (E) = -r sf (E)f 2 ^E)-t nf (E)f n (E), (2) 

where the functions f ncr (E) are the equilibrium distributions of the carriers of a-spin ori- 
entation in region F n , with n = 1 or 2. To make use of these formulae, we observe that 
Boltzmann statistics implies fx± = J2\^~ eV ^ kT ■, and that, as will later be demonstrated in the 
general calculation, the coefficient r 8 f is very small at all energies. We then integrate over 
all energies to obtain the total current in each region, and impose continuity conditions at 

x = and x = d to get 

Js(-d/2) i + i + e~^ kT 

J s (d/2) t+ + t-e- eV / kT [) 

where i± = t n f ± t s f, and the two terms in the sum are population-averaged transmission 

coefficients given by 

f?t nfM (E)e-WdE 
WW) J o °°e~ E / kT dE ' [) 



Together with the standard drift- diffusion theory and other observations noted in Ref. |7[, 
the continuity condition yields the following expressions for the charge current and spin 
currents near the domain wall as functions of voltage and temperature: 



I = Slnh (W 



\2kT 



l + t/ tan h 2 
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~2kT 



— = 2sinh 
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(5) 



D being the diffusion 



where the upper sign holds in F 2 , the lower sign in iq, and Jo 
constant, n< the equilibrium value of the minority spin density, and L s the spin diffusion 
length. Clearly the I — V characteristics of the diode depend critically on the value of the 
i n f/i s f, which will hereafter be referred to as the "key ratio." 

In order to calculate the reflection/transmission probabilities, we must solve the 
Schrodinger equation for the electron wave function in the domain wall 



2m dx 2 



smU{x) cos ti(x) 
cos8(x) — sin6(x) 



+ V{x) 





(6) 



where V(x) = —e£x is the term associated with the electric field £ that is created by 
potential applied across the domain wall. The presence of this term prevents us from finding 
a purely analytical solution, and a numerical solution is therefore computed. Imposing 
the appropriate matching conditions at the domain wall interfaces at x = and d, the 
transmission/reflection probabilities are obtained. When bias produces an electric field 
such that eSd is of the same order as the spin-splitting in this region, the values of these 
probabilities change according to whether the field accelerates or impedes the motion of 
incident electrons through the wall. To assist in observing these effects, we define the 
dimensionless parameters A = ^.^^ which measures the relative size of the domain wall 
barrier, and e = h2 l^md 2 > wn i cn measures the relative strength of the electric field. Values 
for A in the range 0.1 — 0.5 will be considered to describe a thin domain wall, 1 — 5 an 



intermediate size one, and 10 — 50 a thick one. We note for the wall in Ref. 



A ~ 70, 



which is within an order of magnitude of the intermediate size range. Fig. 1(a) shows 



the four coefficients as a function of electron energy for A = 2.25 and zero electric field. 
The essential trends can be easily discerned: at energies less than A, r„/ is approximately 
unity as expected, since the barrier dwarfs the energy of the incident electron. As the 
energy increases, r n f begins to drop and t s f rises at the same rate, since it is now possible 
for the electron to cross the barrier if spin alignment is reversed. At the splitting energy 
threshold, the electron has sufficient energy to traverse the domain wall while maintaining 
spin orientation; t n f increases rapidly while r n f and t s f plummet. We note finally that r s f 
remains approximately zero uniformly over all energies, as previously announced. 

The introduction of an electric field due to current flow has the effect of splitting the 
relevant energy thresholds (Fig. l(b,c)), and the size of the domain wall will determine 
whether this shift is consequential. The minimum energy required for transmission without 
spin-flip is reduced to A — e. The trends follow in a very similar fashion, with the transmis- 
sion/reflection coefficients in the energy range (0,A — e) reaching approximately the same 
values as their zero-field counterparts in the range (0,A), but doing so more rapidly in the 
narrower energy interval, while the coefficients for energy larger than A — e tend to move 
more gradually toward the same limits (£„/ — > 1 and r n f — > as the electron energy E 
grows). Of course electrons of smaller energy can now be transmitted through the reduced 
barrier, thus t n f jumps at this earlier energy threshold, and again at the original barrier 
energy A just slightly. As e exceeds A, the transmission probability is significant at almost 
all non-zero energies; t n f continues to increase uniformly while all other coefficients are sup- 
pressed. This will occur almost immediately for small values of A. For large values of A, 
however, one would have to go to e ^ A in order to have a substantial level of minority spin 
injection: but at this point the resistance of the junction would be too small to support such 
a large electric field. Hence the influence of the electric field is profound for thin domain 
walls and essentially negligible for thick ones. 

These observations account for the main aspects of behavior of the key ratio as a function 
of electric field (see Fig. 2). Physically, values of the key ratio greater than unity signify 
the dominance of minority-spin injection. Again, for A < 0.5, the key ratio is tremendously 
amplified by the electric field, since in this limit t s f goes to zero, and minority-spin injection 
is guaranteed for almost any temperature low enough not to disturb spin-polarization in the 
conductors, but high enough to produce an ample supply of carriers above the exchange 
barrier (this range is typically given by O.lA/k < T < 0.9A/k). The key ratio depends 
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FIG. 1: The reflection and transmission coefficients for an intermediate size domain wall (A = 2.25) 
vs incident electron energy for three values of the electric field : (a) zero field, (b) an electric field 
interaction about half the size of the splitting (e = 1), (c) electric field exceeding the splitting 
(e = 3). The labels of the various curves are shown in (a). 



linearly on temperature for any value of A and e, thus for larger, intermediate barrier sizes, 
there will be a cut-off temperature below which majority spin transmission prevails since 
most of the system's electrons lack sufficient thermal energy to transmit through the wall 
without spin-flip. Our previous observations of effective barrier reduction due to forward 
bias imply that this cut-off temperature will shift downward generally. Indeed, Fig. 2 depicts 
the behavior of the key ratio over a feasible temperature range for A = 5. The zero-field 
curve falls wholly under the minority-spin injection threshold, ~^- = 1, for this barrier size, 
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FIG. 2: (a) Key ratio vs. kT/A for several values of the the electric field parameter, e = 0, 1, ..., 5 
from bottom up at A = 5. The dashed line represents the threshold for minority-spin injection, 
(b) Key ratio vs dimensionless electric field e/A for kT/A = 1 — 10 from bottom up. 



while those for finite values of e exceed it at increasingly lower temperatures: for e = 5 the 
key ratio lies completely above unity. The decaying exponential under the integral in the 
expression for t n f suggests that, for a given value of A, t n f/i s f will rise most rapidly when 
the spin-splitting and the bias voltage have comparable magnitude. This behavior is clearly 
seen in the exponential increase of the key ratio vs electric field in Fig. 2(b), otherwise the 
ratio is approximately linear with voltage for any barrier size. Thus we expect that the 
temperature window of device operation, bounded by the requirements for sufficient carrier 
energy and maintenance of ferromagnetism, will expand downward for intermediate barrier 
sizes, or equivalently, that larger barriers can be accommodated for a fixed temperature 
while still preserving minority-spin injection. 

We are now ready to discuss the behavior of the I — V characteristics, calculated according 
to Eqs. (J5J). Clearly when the key ratio is very large, say > 5, the contribution of the second 
term in the square brackets of Eqs. © is completely negligible. In this case the spin current 
J s reduces to a strictly even function of voltage, the ratio of the spin to the charge current 
Js/Jq (which serves as a measure of spin polarization) is odd-in- volt age, and they are both 
non- linear. Fig. 3 (a-c) shows an example of the behavior of the spin current, the charge 
current, and their ratio, within the first region Fi, The dashed lines show the results 
obtained from Eqs. (J3J) when the value of key ratio is set to the zero-field value. The domain 
is again of intermediate thickness (A = 2.5), but at temperature T = 0.5A/A; the zero-field 
value is clearly quite small and J s has a large odd-in- voltage component. When the voltage 
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FIG. 3: Plots of the spin current (J s , the charge current J q , and their ratio J s /J q vs bias voltage for 
A = 2.5 and two different values of the temperature: (a)-(c) kT/A = 0.5 and (d)-(f) kT/A = 0.2. 
The dashed lines show the results obtained by treating the key ratio -j^- as a constant equal to its 
zero-field value, while the solid line is the result obtained with the voltage-dependent key ratio. 



dependence of the key ratio is included, its rapidly increasing behavior, previously noted, 
leads to a quite different curve, which is shown by the solid line. This is clearly much closer 
to the "ideal" behavior of the spin current, described in Ref. |7]. Similar trends are observed 
for all intermediate values of A, suggesting that, for reasonably sized domain walls, the 
presence of the electric field establishes a new temperature regime in which the key ratio 
is essentially infinite. Yet even outside this regime, it is unlikely for the currents to behave 
as if t n f/t s f were constant. Returning to our example, if the spin splitting energy is kept 
constant but the temperature drops to kT/A = 0.2, then the spin current, as shown in Fig. 
3(d), reaches a negative minimum at a positive voltage determined by the magnitude of A, 
after which it grows again monotonically as a function of voltage. Our study of the key ratio 



suggests that the factors which encourage spin-flip processes in the low-temperature/thick- 
wall regime are eventually overwhelmed by the tendency of the electric field to suppress 
them. While unlikely to contribute to the possible uses of the spin diode, this observation 
does eliminate the possibility of the spin current altering its voltage dependence in a way 
resembling an odd function of voltage. 

In conclusion we have shown that the electric field can assist in maintaining the spin 
polarization of carriers traversing a magnetic domain wall, and consequently the ideal I — V 
characteristics of the spin diode should be more easily attainable than expected. The effect of 
the electric field is also conceptually important since it allows us to better justify the Shockley 
assumptions mentioned earlier in this paper. The reduction in the tunneling probability of 
majority-spin carriers (observed for not too thick domain walls) helps us support assumption 
(1) above. Assumption (3) holds if the spin splitting is large compared to the temperature, 
which can be satisfied more easily now that lower temperatures are acceptable. The increased 
acceleration of the electrons through the domain wall intuitively implies a reduction in the 
time required to traverse its length, which directly supports assumption (4). Finally, we note 
that the foregoing analysis determines the behavior of a junction in terms of the currents 
associated with the four carrier species (spin up and down, n and p). Thus bipolar magnetic 
domain walls and magnetic domain walls between materials with different values of A can be 
analyzed using this approach, simply by using the parameter A associated with the material 
the carriers originate from, and determining e from the band lineup of either the conduction 
bands of the two materials (for electron spin transport) or the valence bands of the two 
materials (for hole spin transport). 
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